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Abstract 

A brief review of checking the factorization hypothesis for matrix element 
of B° — B° mixing within operator product expansion and QCD sum rules 
is given. Both perturbative and power corrections are considered. 
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1. Introduction 

Determination of the exact pattern of CP violation is one of the most important prob- 
lems of modern particle physics. The system of neutral B mesons can provide some useful 
experimental information on the subject. To decipher this information and convert it into 
some knowledge of theoretical parameters of standard model or some extended model 
one needs quite accurate theoretical calculation of corresponding observables within the 
adopted theory. At present an essential obstacle in getting precise theoretical estimates 
for characteristics of CP violation related processes is a necessity of computing hadronic 
matrix elements that is a completely nonperturbative problem. The most popular approx- 
imation for estimating such elements is the factorization, or vacuum saturation, hypoth- 
esis the justification of which is quite unclear. In the present note we very briefly review 
some recent results of checking the validity of the factorization hypothesis with analytical 
methods. We consider three point correlator for computing power corrections violating 
the factorization approximation and two point correlator for computing perturbative ones. 

2. Power corrections: three point correlator 

In order to develop a machinery of operator product expansion and QCD sum rules 1 ' 2 ' 3 
we use here a three point correlator of the form 4,5 

n»v(p,p') = i 2 Jdxdye ipx - ip ' y (0\TJ,(x)O(0) AB=2 J u (y)\0) 

= p^U^p^p'^q 2 ) + ...= Pfl q u U 2 (p 2 ,p' 2 ,q 2 ) + ... (1) 

where q = p' — p, q 2 = 0, J M = d^^^b is an interpolating current for B meson, 

(0\J,\B°(p))=if B p^. 

The operator Oab=2 is chosen with a standard normalization 

A b=2 = &7 M (1 + 75)d&7 M (l + 7 5 )rf 

and the parameter B B is defined by the relation 

(B°\0 AB=2 \B°) = S -f 2 B m 2 B B B . 

Within factorization approximation B B = 1. For higher reliability we take for our analysis 
two invariant functions Ui }2 (p 2 , p' 2 , q 2 ) that appear in the expression for the three point 
correlator Eq. (JJ). The dispersion relation dictates the following representation for these 
functions after using saturation with the low lying resonance state 

IMP ,P ) 1,2=0 = / dsds ~, 5— ~py\ = T~2 2XTT2 2\ + ■ ■ ■ 



[s — p 2 )(s' — p' 2 ) (p 2 — m 2 B )(p' 2 — m 2 B ) 
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For both theoretical and physical parts of sum rules the factorization corresponds to 
the following representation of the amplitudes 

where 

TM = i J dxe^(0\TJ,(x)b(0) lp (l + l5 )dmO), 

and gives the value Bb = 1 for the parameter that describes the ratio of exact matrix 
element to the factorized one. 

Leading contributions of power corrections that are not caught within the factorization 
approximation are given by some specific diagrams with external vacuum fields. The first 
one is given by the gluon condensate 5 

AH? = -^(GG)(pp')(5ri P y(p> 2 ) + e(p 2 )e(p' 2 )), 

An« = AH? - -^(GG)p 2 (r(p 2 )r(p' 2 ) + e{p 2 )e{p' 2 )) + ^ 2 {GG)e{p 2 )g{p' 2 ) 
where 

r(p 2 ) = / dx — - — ■ -, e{p 2 ) = / dx— — 

Jo —p z x + m z Jo —p z x + m z 



g(p 2 ) = — dx(l — 2x)ln(l — xp 2 /m 2 ), 
Jo 



m is the b quark mass, (GG) = (fG^ u G» u ). 

Next contributions are due to condensates of operators with dimension five in mass 
units. They are 



6tt 2 (p 2 — m?) Jo —p' 2 x + m? 

Finally, four quark operators contribute the amount that is not taken into account by the 
factorization approximation for vacuum expectation values 

4 4 2(fo(l+75)^) 2 )"°"- /arf 

An, (p, P ) = An 2 (p, P ) = (p2 _ m2)(p ,_ m2) • 

As in most applications of operator product expansion we do not include operators 
with dimension higher than six in our analysis. Numerical results are obtained after using 
Borel transformation in both p 2 and p' 2 independently and putting M 2 = M' 2 afterwards. 
For both invariant functions there is a fairly wide window of stability with respect to 
change of the Borel parameter. The numerical results are rather stable and reveal only 
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small violation of factorization. Namely, if B B = 1+AB B then for both invariant functions 
and for wide range of parameters (m 2 ,, (qq), ■ ■ ■ ) we get 5 

-AB B = 0-5-0.1. 

This estimate is very conservative. The absolute value of deflection from factorization 
approximation is about —0.05 for the parameter Bb- Main uncertainty is due to poor 
knowledge of numerical value of f B on which there is a strong dependence (to the fourth 
power) . 

3. Perturbative corrections: two point correlator 

Perturbative corrections of order a s that violate the factorization approximation are 
connected with genuine three loop massive diagrams. Their computation with layout for 
sum rules technique based on three point correlator with independent Borel transforma- 
tion with regards to kinematical variables p 2 and p' 2 can not be done at present because 
of technical complexity. Therefore we turn to two point correlator 6,7 and introduce a 
quantity 8 

T(x) = (0\TO AB=2 (x)O AB=2 (0)\0). 

The leading term of a s expansion for the above correlator has an expression in the con- 
figuration space that reads 

T (x) = 2N 2 (l + 16S'(x,m)S(-x,0)S'(x,m)S(-x,0) 

= 2 (l + ^) tr[S(x, m)S(-x, 0)}tr[S(x, m)S(-x, 0)] = 2 (l + i-) U 5 (x)U 5 (x) (2) 

where S(x, m) is the free fermion propagator and N c stands for the number of quark 
colors. The prime means taking only the part of the propagator that is proportional to a 
7 matrix. The function n 5 (a;) = (0\Tj 5 (x)j 5 (0)\0) is the two point correlator associated 
to the current j 5 = bi^d. Thus one observes a complete factorization in this order. 
Eq. (2) can be rewritten in the form 

T (x) = 2{l + j^jU, v (x)W^x), (3) 

where W v (x) = (0\Tj£(x)f L (0)\0) and f L = b L ^d L . The Lorentz decomposition in 
x-space reads 

W(x) = (-d»d v + g" u d 2 )U T {x 2 ) - d^d"U L {x 2 ) 

that again demonstrates an explicit factorization in the configuration space to leading 
order in a s . 
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The dispersion representation in x-space for any two point correlator IL(x) (j = 
T, L, 5) has the form 

ill j (x 2 )= rj(s)D(x,s)ds (4) 



where D(x, s) is a free boson propagator with the "mass" y/s. The spectral functions r,- 
read to leading order in a s 

^M^ 2 ' 1 -^ 4°»W = ^(1-^(2 + ,), r fM = ^^ 

where z = m 2 /s, m is the b quark mass. 

The spectral function p(s) of the full correlator T(x) is defined in the same way as in 
Eq. (4). To first order in a s it can be expressed in terms of the spectral functions rj(s) 
associated to the two-line correlators in the form 

p(s) = J ri(si)r 2 (s 2 )$(s; s u s 2 )ds 1 ds 2 

where 

si, s 2 ) = 16?r 2g \l (~ s ~ Sl ~ g2 ) 2 ~ 4glg2 

is the two-body phase space factor. 

To leading order in 1/N C , one can write the correlator T(x) as a product of two two- 
line correlators. It is worthwhile to notice that this decomposition is gauge invariant and 
finite, i.e. it does not require any renormalization. 

The full spectral density p(s) has been computed numerically to the first order in 
a s 8 with a heavy use of known results for two loop massive diagrams obtained earlier 
(e.g. 9 ) and the program of symbolic computation REDUCE. We analyze our results con- 
centrating on presentation of the entire spectral density p(s) as a sum of factorizable and 
nonfactorizable pieces 

p(s) = p (s) (1 + Ap f (s) + A Pnf (s)) . 

Nonfactorizable part of the spectral density is given by one gluon exchange diagrams of 
a two-line correlator. 

Within sum rules approach one works with moments 

rs t h 

Mi(s th ) = / p(s)s- l ds 

J4m 2 b 

that are decomposed as 

Mi(sth) = M?(s th ) (l + AM/ (s th ) + AM? f (s th j) 
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according to the decomposition of the spectral density. 

A set of input parameters for numerical estimates is =175 MeV, m b = 4.6 GeV. 
The moments of the factorizable spectral density are almost independent of the power i 
of the weight function s~\ The nonfactorizable correction does not exceed a 15% level 
with respect to the full factorized spectral density. Corrections of order a s (factoriz- 
able+nonfactorizable) are large for both the spectral density itself and its moments. De- 
pending on the energy s they can reach a magnitude of 100% with respect to the leading 
term. The nonfactorizable corrections measured in terms of the fully factorized (lowest 
order + a s terms) spectral density are moderate. All the factorizable corrections to the 
correlator, however, can be absorbed into the calculation of the decay constant fs from 
two point correlator with two quark lines, in such a way that the relevant corrections to 
the Bb parameter are only due to nonfactorizable ones. Results for spectral density itself 
and for its moments for different values of integration regions are collected in Tables 1,2 
taken from ref. 8 . 
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A Pf 


Ap nf 


A Pnf /(l + A Pf ) 


5.5 


1.03 


0.02 


0.01 


6.0 


0.95 


0.21 


0.11 


6.5 


0.89 


0.29 


0.15 


7.0 


0.84 


0.32 


0.17 



Table 1: Normalized spectral densities 



% 




AM f 


AM nf 


AM n/ /(l + AMf) 





5.5 


1.07 


-0.16 


-0.08 




6.0 


0.99 


0.11 


0.06 




6.5 


0.93 


0.23 


0.12 




7.0 


0.88 


0.29 


0.15 


5 


5.5 


1.08 


-0.21 


-0.10 




6.0 


1.00 


0.08 


0.04 




6.5 


0.94 


0.20 


0.11 




7.0 


0.89 


0.27 


0.14 


10 


5.5 


1.09 


-0.27 


-0.13 




6.0 


1.01 


0.02 


0.01 




6.5 


0.96 


0.16 


0.08 




7.0 


0.91 


0.23 


0.12 



Table 2: Normalized moments of spectral densities 
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